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Devices for nano- and molecular size electronics are currently a focus of 
research aimed at an efficient current rectification and switching. Current 
switching due to conformational changes in the molecules is slow, on the order 
of a few kHz. Fast switching (^1 THz) may be achieved, at least in princi- 
ple, in a degenerate molecular quantum dot with strong coupling of electrons 
with vibrational excitations. We show that the mean-field approach fails to 
properly describe intrinsic molecular switching and present an exact solution 
to the problem. 

1 Introduction 

For many applications one needs an intrinsic molecular "switch" , i.e. a bistable 
voltage-addressable molecular system with very different resistances in the two 
states that can be accessed very quickly pQ. There is a trade-off between the 
stability of a molecular state and the ability to switch the molecule between 
two states with an external perturbation (we discuss an electric field, switch- 
ing involving absorbed photons is impractical at a nanoscale). Indeed, the ap- 
plied electric field, on the order of a typical breakdown field Ei, < lO^V/cm, 
is much smaller than a typical atomic field ~ lO^V/cm, characteristic of the 
energy barriers. Small barrier would be a subject for sporadic thermal switch- 
ing, whereas a larger barrier ~ 1 — 2eV would be impossible to overcome 
with the applied field. One may only change the relative energy of the min- 
ima by external field and, therefore, redistribute the molecules statistically 
slightly inequivalcntly between the two states. An intrinsic disadvantage of 
the conformational mechanism 2 , involving motion of ionic group, exceeding 
the electron mass by many orders of magnitude, is a slow switching speed 
(~kHz). In case of supramolecular complexes like rotaxanes and catenanes |3] 
there are two entangled parts which can change mutual positions as a result 
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Fig. 1. Schematic of tiie molecular quantum dot with central conjugated unit sep- 
arated from the electrodes by wide-band insulating molecular groups. First electron 
tunnels into the dot and occupies an empty (degenerate) state there. If the inter- 
action between the first and second incoming electron is repulsive, (7 > 0, then the 
dot will be in a Coulomb blockade regime (a) . If the electrons on the dot effectively 
attract each other, 17 < 0, the system will show current hysteresis (b). 



of redox reactions (in solution). Thus, for the rotaxane-based memory devices 
a slow switching speed of ^ 10~^ seconds was reported. 

We have, therefore, explored a possibiUty for a fast molecular switching 
where switching is due to strong correlation effects on the molecule itself, 
so-called molecular quantum dot (MQD). The molecular quantum dot con- 
sists of a central conjugated unit (containing half-occupied, and, therefore, 
extended tt— orbitals), Fig. ^ Frequently, those are formed from the p-states 
on carbon atoms, which are not saturated (i.e. they do not share electrons 
with other atoms forming strong a— bonds, with typical bonding-antibonding 
energy difference about IRy). Since the tt— orbitals are half-occupied, they 
form the HOMO-LUMO states. The size of the HOMO-LUMO gap is then 
directly related to the size of the conjugated region d. Fig. ^ by a standard 
estimate -Ehomo-lumo ^ h^/md'^ ~ 2 — 5 eV. It is worth noting that in the 

conjugated linear polymers like polyacetylene (— C=C)ri the spread of the 
TT— electron would be d = oo and the expected ii^noMO-LUMO = 0- However, 
such a one-dimensional metal is impossible, Peierls distortion (C=C bond 
length dimerization) sets in and opens up a gap of about ~ 1.5eV at the 
Fermi level In a molecular quantum dot the central conjugated part 

is separated from electrodes by insulating groups with saturated cr— bonds, 
like e.g. the alkane chains, Fig. 3. Now, there are two main possibilities for 
carrier transport through the MQD. If the length of at least one of the insu- 
lating groups ii(2) is not very large (a conductance Gi(2) is not much smaller 
than the conductance quantum Go = 2e^//i), then the transport through the 
MQD will proceed by resonant tunneling processes. If, on the other hand. 
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both groups are such that the tunnel conductance Gi(2) ^ Go, the charge on 
the dot will be quantized. Then we will have another two possibilities: (i) the 
interaction of the extra carriers on the dot is repulsive U > 0, and we have a 
Coulomb blockade 0, or (ii) the effective interaction is attractive, U < 0, then 
we would obtain the current hysteresis and switching [S] (see below). Coulomb 
blockade in molecular quantum dots has been demonstrated in Refs. 0. In 
these works, and in Ref. the three-terminal active molecular devices have 
been fabricated and successfully tested. 

Much faster switching compared to the conformational one may be caused 
by coupling to the vibrational degrees of freedom, if the vibron-mediated 
attraction between two carriers on the molecule is stronger than their direct 
Coulomb repulsion, Fig.nj). The attractive energy (i.e. a negative "Hubbard" 
U) is the difference of two large interactions, the Coulomb repulsion and the 
phonon mediated attraction, on the order of leV each, hence \U\ ^ O.leV. 

2 Failure of mean field model of polaron molecular 
switching 

Although the correlated electron transport through mesoscopic systems with 
repulsive electron-electron interactions received considerable attention in the 
past, and continues to be the focus of current studies, much less has been 
known about a role of electron-phonon correlations in "molecular quantum 
dots" (MQD). Some while ago we have proposed a negative— f7 Hubbard 
model of a d-fold degenerate quantum dot ^Jj and a polaron model of reso- 
nant tunneling through a molecule with degenerate level [H| . We found that the 
attractive electron correlations caused by any interaction within the molecule 
could lead to a molecular switching effect where I-V characteristics have two 
branches with high and low current at the same bias voltage. This prediction 
has been confirmed and extended further in our theory of correlated transport 
through degenerate MQDs with a full account of both the Coulomb repulsion 
and realistic electron-phonon (e-ph) interactions. We have shown that while 
the phonon side-bands significantly modify the shape of hysteretic I-V curves 
in comparison with the negative- 1/ Hubbard model, switching remains robust. 
It shows up when the effective interaction of polarons is attractive and the 
state of the dot is multiply degenerate, d > 2. 

Nevertheless, later on Galperin et al. argued, without discussing the 
discrepancies with the prior work, that even a non-degenerate electronic level 
{d = 1) coupled to a single vibrational mode produces a hysteretic TV curve, 
a current switching, and a negative differential resistance. Here we explicitly 
calculate TV curves of the nondegenerate {d = 1) and two-fold degenerate 
{d = 2) MQDs to show that these findings are artefacts of the mean-field 
approximation used in Ref. ^21 that neglects the Fermi-Dirac statistics of 
electrons. 
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Fig. 2. Two localized electrons at sites 1 and 2 shift the equilibrium position of the 
ion at site 3. As a result, the two electrons attract each other. 



We start with a simple model that illustrates an absence of switching 
in molecular quantum dot, which has non-degenerate {d — 1) or double- 
degenerate (d = 2) level First, we shall illustrate the failure of the mean- 
field approximation on a simplest model of a single atomic level coupled with a 
single one-dimensional oscillator with a displacement a;, described by a simple 
Hamiltonian, 

H = e.n + f.n-——^+—. (1) 

Here M and k are the oscillator mass and the spring constant, / is the inter- 
action force, and h — c = ks — This Hamiltonian is readily diagonalized 
with the exact displacement transformation of the vibration coordinate x, 

x^y- nf/k, (2) 

to the transformed Hamiltonian without electron-phonon coupling, 

H = en H , (3) 

2M dy^ 2 ' ^ ' 

e = eo - Ep, (4) 

where we used fi^ — h because of the Fermi-Dirac statistics. It describes a 
small polaron at the atomic level Eq shifted down by the polaron level shift 
Ep — f'^ /2k, and entirely decoupled from ion vibrations. The ion vibrates near 
a new equilibrium position, shifted by f/k, with the "old" frequency (fc/M)^/^. 
As a result of the local ion deformation, the total energy of the whole system 
decreases by Ep since a decrease of the electron energy by —2Ep overruns 
an increase of the deformation energy Ep. The major error of the mean- field 
approximation of Ref. ,12, originates in illegitimate replacement of the square 
of the occupation number operator n = cJcq by its "mean- field" expression 
= nofi which contains the average population of a single molecular level, ng , 
in disagreement with the exact identity, — h. This leads to a spurious self- 
interaction of a single polaron with itself [i.e. the term e = e^ — n^Ep instead 
of Eq. (0}], and a resulting non-existent nonlinearity in the rate equation. 
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Lattice deformation also strongly affects the interaction between electrons. 
When a short-range deformation potential and molecular e-ph interactions are 
taken into account together with the long-range Frohlich interaction, they can 
overcome the Coulomb repulsion. The resulting interaction becomes attrac- 
tive at a short distance comparable to a lattice constant. The origin of the 
attractive force between two small polarons can be readily understood from a 
similar Holstein-like toy model as above ^14, , but with two electrons on neigh- 
boring sites 1,2 interacting with an ion 3 between them, Fig. 12 For generality, 
we now assume that the ion is a three-dimensional oscillator described by a 
displacement vector u, rather than by a single-component displacement x as 
in Eq.(l). 

The vibration part of the Hamiltonian in the model is 
„ _ I fcu2 

Electron potential energies due to the Coulomb interaction with the ion are 
about 

^1,2 = l^o(l-u-ei,2/a), (6) 

where ei 2 are units vectors connecting sites 1,2 and site 3, respectively. 
Hence, the Hamiltonian of the model is given by 

H = Ea{ni + n2) + u • (fi^i + hn^) " TTTFIT^ + (7) 

where fi.2 = Ze^Gx^-ij a} is the Coulomb force, and n\,i are occupation number 
operators at every site. This Hamiltonian is also readily diagonalized by the 
same displacement transformation of the vibronic coordinate u as above, 

u = v-(firii+f2n2)/fc. (8) 

The transformed Hamiltonian has no electron-phonon coupling, 

19^ fcv^ 

H ^ {eo- Ep){hi + na) + Vphnin2 " + — ' (9) 

and it describes two small polarons at their atomic levels shifted by the polaron 
level shift Ep — fi2/'2k, which are entirely decoupled from ion vibrations. As 
a result, the lattice deformation caused by two electrons leads to an effec- 
tive interaction between them, Vpi-^ which should be added to their Coulomb 
repulsion, Vc, 

Vph = -fi • h/k. (10) 

When Vph is negative and larger by magnitude than the positive Vc, the re- 
sulting interaction becomes attractive. That is Vph rather than i?p, which is 
responsible for the hysteretic behavior of MQDs, as discussed below. 
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3 Exact solution of polaron switching 

The procedure, which fuUy accounts for all correlations in MQD is as follows, 
see Ref. [HI- The molecular Hamiltonian includes the Coulomb repulsion, f/*^, 
and the electron- vibron interaction as 



F = y^e?i +-y^ U^n'n flu' 

+ f'.M^f^idg + H.c.) + '^M^i + 1/2). (11) 

Here dq annihilates phonons, Uq is the phonon (vibron) frequency, and 
are the electron- vibron coupling constant {q enumerates the vibron modes). 
This Hamiltonian conserves the occupation numbers of molecular states h^. 
One can apply the canonical unitary transformation e"^, with 

S" = - ^ 71^ {jf,qdq - H.c.) 

integrating phonons out. The electron and phonon operators are transformed 
as 



(12) 



and 

dq ^ dq-J2f^t^%q^ (13) 

respectively. This Lang-Firsov transformation shifts ions to new equilibrium 
positions with no effect on the phonon frequencies. The diagonalization is 
exact: 

H = + II^«(4^9 + 1/2) + \Y. '^mm'^mV. (14) 

where 

U^.^.' = C/^^' - 2 ^ %ql^.'qUJq, (15) 

q 

is the renormalized interaction of polarons comprising their interaction via 
molecular deformations (vibrons) and the original Coulomb repulsion, C/^/ . 
The molecular energy levels are shifted by the polaron level-shift due to the 
deformation created by the polaron, 



^.-El^M^l''^?- (16) 



9 



If we assume that the coupling to the leads is weak, so that the level width 
J" <C \U\, we can find the current from 
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i{v)^er d^[/iH-/2H]pH 

J —oc 



(17) 
(18) 



where |/Lt) is a complete set of one-particle molecular states, fi(2){^) 



expi^±^i^ + l 



the Fermi function. Here p{lo) is the molecular DOS, 



( 

G^{lo) is the Fourier transform of the Green's function G^{t) = —i9{t) (|cp(t), cj,}) 
{• • • , • • • } is the anticommutator, c^{t) = e'^*c^e~*^*, 9{t) — 1 for t > and 
zero otherwise. We calculate p{u!) exactly for the Hamiltonian p4|) . which 
includes both the Coulomb and electron-vibron interactions. 
The retarded GF becomes 



G^{t)^^^0{t)[{c,{t)cl){X,{t)Xl) 
+ {clc,,{t)){XlX,{t))]. 

The phonon correlator is simply 



(19) 



{X,{t)Xl)^expJ2 



I 7m 



sinh 



cos [ ujt + i 



— cosh 

2 



(20) 



where the inverse temperature (3 = 1/T, and (^Xj^X^(t)'^ = (^X^{t)Xj^^ ' . The 
remaining GFs (^Ci^(t)cjj^, are found from the equations of motion exactly. 
For the simplest case of a coupling to a single mode with the characteristic 
frequency wq, 7^ =7 and f/^^/ = U one obtains |H| 



G»=Z^a(n)^/KO 



r=0 



e 2 



1 



n 



n 



+ (1 -<5;o)e 
1 



uj — rU — Iluq + iS uj — rU + Iluq + i5 



n 



n 



where 



UJ — rU + lujQ + i5 Lo — rU — Iljq + iS 



Z = exp ^-X! I^J^coth^^ 



is the familiar polaron narrowing factor, the degeneracy factor 

{d-iy. 



Cr{n) 



rl{d - 1 - r)! 



(21) 



(22) 



(23) 
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C = I7P/ sinh^^^, /; (^) the modified Bessel function, and 6i^. the Kroneker 
symbol. 

Then using Eq.(18) the exact spectral function for a d — fold degenerate 
MQD (i.e. the density of molecular states, DOS) is found as 



d-l 



r=0 



X e'^""'/^ [(1 - n)d{LU -rU - luJa) + nS{uj - rU + Iloq) 
+ (1 - Sio)e-^'^°^/^[nS{uj - rU - Iujq) 
+ (1 - n)5{Lo -rU + IijJo)] 



(24) 



The important feature of DOS, Eq. (|24|l . is its nonlinear dependence on 
the average electronic population n — l^c^^c^^ , which leads to the switching, 
hysteresis, and other nonlinear effects in I-V characteristics for d > 2. It ap- 
pears due to correlations between different electronic states via the correlation 
coefficients Cr{n). There is no nonlinearity if the dot is nondegenerate, d—\, 
since Co(n) = 1. In this simple case the DOS, Eq. H24|l . is a linear function of 
the average population that can be found as a textbook example of an exactly 
solvable problems 16 . 

In the present case of MQD weakly coupled with leads, one can apply 
the Fermi-Dirac golden rule to obtain an equation for n. Equating incoming 
and outgoing numbers of electrons in MQD per unit time we obtain the self- 
consistent equation for the level occupation n as 



= n 



(l-n) dLu{rifi{i0)+r2f2{L0)}pii0) 
J — 00 

{ A [1 - /l (C^)] + 1^2 [1 - /2 (C^)] } p(cc') , (25) 



— 00 



where A(2) are the transition rates from left (right) leads to MQD, and p{llj) 
is found from Eq.(24). For d = 1,2 the kinetic equation for n has only one 
physical root, and the switching is absent. Switching appears for d > 3, when 
the kinetic equation becomes non-linear. Taking into account that J_ Pi'^) = 
d, Eq. H25|l for the symmetric leads, A = /2, reduces to 

2nd = J duop iu) (/i + /2) , (26) 

which automatically satisfies < n < 1. Explicitly, the self-consistent equation 
for the occupation number is 

n = - Zr{n)[nar + (1 — n)br], (27) 

r=0 
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where 

OO y 

a+ = Z ^ (0 I [/i {rU - lu^) + h {rU - lu^)] 

1=0 ^ 

+ (1 - Sio)e-^ [h{rU + Icuo) + f2{rU + lujo)]^ , (28) 

OO / 

1=0 ^ 

+ (1 - 5iQ)e-^ [h{rU - Icuo) + f2{rU - lu;o)]j . (29) 
The current is expressed as 

= ^ = E Zrin)[na- + (1 - n)b;], (30) 
-'0 



where 



1=0 



r=0 



a- = Z ^ Ii (0 ( [MrU - Iujq) - f2{rU - loJo)] 



+ (1 - Sio)e-— [f,{rU + Itoo) - f2{rU + luo)]j , (31) 

OC ^ 

b- =ZY, h (0 ( [/i {rU + luo) - /2(r[/ + luj^)] 
1=0 ^ 

(1 - 5io)e-^ [h{rU - Zc^o) - f2{rU - luj^)]] , (32) 



and Iq =^ edr. 



4 Absence of switching of single- or double-degenerate 
MQD 

If the transition rates from electrodes to MQD are small, F <^ luq, the rate 
equation for n and the current, I{V) are readily obtained by using the exact 
molecular DOS, Eq. (|24|l and the Fermi-Dirac Golden rule. In particular, for 
the nondegenerate MQD and T = OK the result is 

n= '$ -, (33) 

and 
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Fig. 3. Current- voltage characteristic of the nondegenerate (d = 1) MQD at T = 0, 
wo/zi = 0.2, and 7^ = 11/13. There is the phonon ladder in I-V , but no hysteresis. 
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Fig. 4. Current-voltage characteristic of two-fold degenerate MQDs {d = 2) does 
not show hysteretic behavior. The parameters are the same as in Fig. El Lar ger 
number of elementary processes for conductance compared to the nondegenerate 
case of d = 1 generates more steps in the phonon ladder in comparison with Fig. El 



■ ^ +a„b+ -a+b^ 
2 + b+-a+ 

The general expressions for the coefficients Eqs. (28,29) and Eqs. (31,32) at 
arbitrary temperatures in Ref. |S] are simplified in low-temperature limit as 

00 I |2Z 

4 = ^[^(^^0 - ^ + eV/2) 

1=0 

±e{lujo-A-eV/2)l (35) 

b^=Zj2^-^m-l^o-^ + eV/2) 
1=0 

±e{-lujo ~ A~ eV/2)], (36) 
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where A is the position of the MQD level with respect to the Fermi level at 
V = 0, and 0{x) = 1 if a; > and zero otherwise. The current is single valued, 
Fig. 131 with the familiar steps due to phonon-side bands. 

On the contrary, the mean-field approximation (MFA) leads to the op- 
posite conclusion. Galperin et al. ^21 have replaced the occupation number 
operator fi in the e-ph interaction by the average population no [Eq. (2) of 
Ref. 521] and found the average steady-state vibronic displacement {d + d^) 
proportional to no (this is an explicit neglect of all quantum fluctuations on 
the dot accounted for in the exact solution). Then, replacing the displacement 
operator d + d'' in the bare Hamiltonian, Eq. (11), by its average, Ref. |12| . 
they obtained a new molecular level, Eq — Eq — 2ereorgno shifted linearly 
with the average population of the level. This is in stark disagreement with 
the conventional constant polaronic level shift, Eq. H4I16|I (Sreorg is ItI^^o in 
our notations). The MFA spectral function turned out to be highly nonlin- 
ear as a function of the population, e.g. for the weak-coupling with the leads 
p{uj) — 5{uj — eo — 2ereorgno), see Eq. (17) in Ref. J^l- As a result, the au- 
thors of Ref. [T2j have found multiple solutions for the steady-state population, 
Eq. (15) and Fig. 1, and switching. Fig. 4 of Ref. which actually do not 
exist being an artefact of the approximation. 

In the case of a double-degenerate MQD, d = 2, there are two terms, 
which contribute to the sum over r, with Co(n) = 1 — n and Ci(n) = n. 
The rate equation becomes a quadratic one |H] ■ Nevertheless there is only one 
physical root for any temperature and voltage, and the current is also single- 
valued. The double-degenerate level provides more elementary processes for 
conductance reflected in larger number of steps on phonon ladder compared 
to d = 2 case, Fig. U] 

Note that the mean-field solution by Galperin et al. 12 applies at any 
ratio r/ujQ^ including the limit of interest to us, <^ wo- where their transition 
between the states with no = and 1 only sharpens, but none of the results 
change. Therefore, MFA predicts a current bistability in the system where it 
does not exist at d = 1. Ref. ^21 plots the results for F > wq, « 0.1 — 0.3 
eV, which corresponds to molecular bridges with a resistance of about a few 
lOOKi?. Such model "molecules" are rather "metallic" in their conductance 
and could hardly show any bistability at all because carriers do not have 
time to interact with vibrons on the molecule. Indeed, taking into account the 
coupling with the leads beyond the second order and the coupling between the 
molecular and bath phonons could hardly provide any non-linearity because 
these couplings do not depend on the electron population. This rather obvious 
conclusion for molecules strongly coupled to the electrodes can be reached in 
many ways, see e.g. a derivation in Refs. J3 ^j. While Refs. ^| E| do 
talk about telegraph current noise in the model, there is no hysteresis in the 
adiabatic regime, F ^ uq either. This result certainly has nothing to do 
with our mechanism of switching j2| that applies to molecular quantum dots 
{F 4C Wo) with d > 2. Such regime has not been studied in Refs. |,17, 18, iHI, 
which have applied the adiabatic approximation, as being "too challenging 
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problem" . Nevertheless, Mitra et al. |19j have misrepresented our formalism 
|S] claiming that it "lacks of renormanlization of the dot-lead coupling" (due to 
electron- vibron interaction), or "treats it in an average manner". In fact, the 
formalism [5! is exact, fully taking into account the polaronic renormalization, 
phonon-side bands and polaron-polaron correlations in the exact molecular 
DOS, Eq. lO. 

As a matter of fact, most of the molecules are very resistive, so the actual 
molecular quantum dots are in the regime we study, see Ref.|^. For exam- 
ple, the resistance of fully conjugated three-phenyl ring Tour-Reed molecules 
chemically bonded to metallic Au electrodes exceeds 1GJ7. Therefore, most 
of the molecules of interest to us are in the regime that we discussed, not that 
of Refs-HailHl. 

5 Nonlinear rate equation and switching 

The switching appears only for d > 2. For example, for d = 4 the rate equation 
(27) is of the fourth power in n, 

2n = (1 - nf[na+ + (1 - n)&+] 

-|-3n(l - nf[naf + (1 ~ n)bt] 
+3n'^{l - n)[na+ + (1 - n)b+] 
+n^[na+ + (l - n)b+]. (37) 

Different from the non-degenerate or double-degenerate MQD, the rate 
equation for d = 4 has two stable physical roots in a certain voltage range and 
the current-voltage characteristics show a hysteretic behavior. Our numerical 
results |H] for uiq = 0.2 (in units of A, as all the energies in the problem), 
C/*^ = 0, and for the coupling constant, 7^ — 11/13 are shown in Fig. 5. 
This case formally corresponds to a negative Hubbard U — —2j'^ljo « —0.4 
(we selected those values of 7^ to avoid accidental commensurability of the 
correlated levels separated by U and the phonon side-bands) . The threshold for 
the onset of bistability appears at a voltage bias eV/2A — 0.86 for 7^ = 11/13 
and wq = 0.2. The steps on the I-V curve. Fig. 5, are generated by the phonon 
side-bands originating from correlated levels in the dot with the energies A, 
A + U, A + {d — 1)U. Since cuo is not generally commensurate with U, we 
obtain quite irregular picture of the steps in I-V curves. The bistability region 
shrinks down with temperature. 

Note that switching required a degenerate MQD {d > 2) and the weak 
coupling to the electrodes, F <^ luq. Different from the non-degenerate dot, 
the rate equation for a multi-degenerate dot, d > 2, weakly coupled to the 
leads has multiple physical roots in a certain voltage range and a hysteretic 
behavior due to correlations between different electronic states of MQD. 
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Fig. 5. The bistable I-V curves for tunneling through molecular quantum dot with 
the electron-vibron coupling constant 7^ = 11/13 and coo/^ ~ 0.2. The up arrows 
show that the current piclcs up at some voltage when it is biased, and then drops 
at lower voltage when the bias is being reduced. The bias dependence of current 
basically repeats the shape of the level occupation n (right column). Steps on the 
curve correspond to the changing population of the phonon side-bands 



6 Summary 

We have calculated the I-V characteristics of a nondegenerate (d = 1) , two-fold 
degenerate {d — 2) molecular quantum dots showing no hysteretic behavior 
of current, and concluded that mean field approximation ^21 leads to a non- 
existent swtching in a model that was solved exactly in Ref. W- Different 
from the non-degenerate and two-fold degenerate dots, the rate equation for a 
multi-degenerate dot, d> 2, weakly coupled to the leads, has multiple physical 
roots in a certain voltage range showing hysteretic behavior due to correlations 
between different electronic states of MQD 8 . Pair tunneling is also allowed 
in out model, though it should only result in tiny peaks on the background 
of the main current contributed by single polaron tunneling. Our conclusions 
are important for searching for the current-controlled polaronic molecular-size 
switches. Incidentally, Cgo molecules have the degeneracy d = 6 of the lowest 
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unoccupied level, which makes them one of the most promising candidate 
systems, if the weak-coupling with leads is secured. 

We thank the participants of the ESF workshop "Mott's Physics in 
Nanowires and Quantum Dots" (Cambridge, UK, 31 July-2 August, 2006) for 
many discussions, and greatly appreciate the financial support of the Euro- 
pean Science Foundation (ESF) and EPSRC (UK) (grants no EP/C518365/1 
and no EP/D07777X/1). 



References 

1. Z.J. Donhauser, B.A. Mantooth, K.F. Kelly, L.A. Bumm, J.D. Monnell, J.J. 
Stapleton, D.W. Price, Jr., A. M. Rawlett, D. L. AUara, J. M. Tour, and P. S. 
Weiss, Science 292, 2303 (2001); Z. J. Donhauser, B.A. Mantooth, T.P. Peari, 
K.F. Kelly, S.U. Nanayakkara, and P.S. Weiss, Jpn. J. Appl. Phys. 41, 4871 

(2002) . 

2. C. Li, D. Zhang, X. Liu, S. Han, T. Tang, C. Zhou, W. Fan, J. Koehne, J. Han, 
M. Meyyappan, A.M. Rawlett, D.W. Price, and J.M. Tour, Appl. Phys. Lett. 
82, 645 (2003); M.A. Reed, J. Chen, A.M. Rawlett, D.W. Price, and J. M. Tour, 
Appl. Phys. Lett. 78, 3735 (2001). 

3. CP. Collier, E. W. Wong, M. Belohradsky, F. M. Raymo, J. F. Stoddart, P. J. 
Kuekes, R. S. WiUiams, and J. R. Heath, Science 285, 391 (1999); CP. Collier, 
G. Mattersteig; E.W. Wong, Y. Luo, K. Beverly, J. Sampaio, P.M. Raymo, J.F. 
Stoddart, J.R. Heath, Science 289, 1172 (2000). 

4. Y. Chen, D.A.A. Ohlberg, X. Li, D.R. Stewart, and R.S. Williams, J.O. Jeppe- 
sen, K.A. Nielsen, and J.F. Stoddart, D.L. Olynick, and E. Anderson, Appl. 
Phys. Lett. 82, 1610 (2003); Y. Chen, G.Y. Jung, D.A A Ohlberg, X. Li, D.R 
Stewart, J.O Jeppesen, K.A Nielsen, J F. Stoddart, and R S. Williams, Nan- 
otechnology 14 462 (2003). 

5. D. R. Stewart, D. A. A. Ohlberg, P. A. Beck, Y. Chen, R. S. Williams, J. O. 
Jeppesen, K. A. Nielsen, and J. F. Stoddart, Nanoletters 4, 133 (2004). 

6. Molecular Switches, edited by Ben L. Feringa (Wiley- VCH, 2001). 

7. D.V. Averin and K.K. Likharev, in: Mesocopic Phenomena in Solids, edited by 
B.L. Altshuler et al. (North-Holland, Amsterdam, 1991). 

8. A.S. Alexandrov and A.M. Bratkovsky, Phys. Rev. B 67, 235312 (2003). 

9. H. Park, J. Park, A.K.L. Lim, E.H. Anderson, A.P. Alivisatos, and P.L. McEuen, 
Nature 407, 57 (2000); J. Park, A.N. Pasupathy, J.L Goldsmith, C. Chang, Y. 
Yaish, J.R. Petta, M. Rinkoski, J.P. Sethna, H.D. Abruha, P.L. McEuen, and 
D.C Ralph, ibid. 417, 722 (2002); W. Liang, M.P. Shores, M. Bockrath, J.R. 
Long, and H. Park, ibid. 417, 725 (2002). 

10. N.B. Zhitenev, H. Meng, and Z. Bao, Phys. Rev. Lett. 88, 226801 (2002). 

11. A.S. Alexandrov, A.M. Bratkovsky, and R.S. WiUiams, Phys. Rev. B 67, 075301 

(2003) . 

12. M. Galperin, M.A. Ratner, and A. Nitzan, Nano Lett. 5, 125 (2005). 

13. A.S. Alexandrov and A.M. Bratkovsky, cond- mat/0606366 

14. A.S. Alexandrov and N.F. Mott, Polarons and Bipolarons (World Scientific, 
Singapore, 1996). 

15. Y. Meir and N.S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 



Polaron switching in molecular devices 15 



16. G.D. Mahan, Many-Particle Physics, 2nd ed. (Plenum, New York, 1993). 

17. A. Mitra, I. Aleiner, and A. Millis, Phys. Rev. Lett. 94, 076404 (2006). 

18. D. Mozyrsky, M. B. Hastings, and I. Martin, Phys. Rev. B 73, 035104 (2006). 

19. A. Mitra, I. Aleiner, and A. Milhs, Phys. Rev. B 69, 245302 (2004). 

20. H. Park, J. Park, A.K.L. Lim, E.H. Anderson, A.P. Alivisatos, and P.L. McEuen, 
Nature 407, 57 (2000); J. Park, A.N. Pasupathy, J.I. Goldsmith, C. Chang, 
Y.Yaish, J.R. Retta, M. Rinkoski, J. P. Sethna, H.D. Abruna, P.L. McEuen, and 
D.C. Ralph, Nature (London) 417, 722 (2000); W. Liang et al, Nature (London) 
417, 725 (2002). 



